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Abstract 

A generalized quantum kinetic equation (RKE) of the Kadanoff-Baym type is ob- 
tained on the basis of the Heisenberg equations of motion where the time evolution and 
space translation are separated from each other by means of the covariant method. The 
same approach is used also for a covariant modification of the real-time Green's functions 
method based on the Wigner representation. The suggested approach does not contain 
arbitrariness' elements and uncertainties which often arise from derivation of RKE on the 
basis of the motion equations of the Kadanoff-Baym type for the correlation functions in 
the case of systems with inner degrees of freedom. 

Possibilities of the proposed method are demonstrated by examples of derivation of 
RKE of the Vlasov type and collision integrals of the Boltzmann-Uehling-Uhlenbeck 
(BUU) type in the frame of the crw-version of quantum hadrodynamics, for the sim- 
plest case of spin saturated nuclear matter without antinuclear component. Here, the 
quasiparticle approximation in a covariant performance is used. A generalization of the 
method for the description of strong non-equilibrium states based on the non-equilibrium 
statistical operator is then proposed as well. 



1 Introduction 

The recent growing interest in the relativistic nuclear physics, especially at intermediate 
energies, has stimulated the development of the theory of relativistic kinetic equations 
(RKE) as well. In particular, the main problems of dynamical RKE derivation were 
discussed within the framework of different approaches with well grounded traditions in 
the non-relativistic case, namely the method of the real-time (contour) Green's functions 
the BBGKY method (method of many particles correlation functions) [0, the 
method of the non-equilibrium statistical operator |9|, and several other techniques 



However, in some approaches, the transition to the relativistic region it is not obvious 
and can lead to misleading results. This situation was analyzed in Refs.[16|-[1?] on a 



rather simple quantum field model which describes interaction of the fermion and boson 
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subsystems in the mean field approximation. In these works, it was shown that the 
equations of motion for a two-point correlation function lead to true RKE of the Vlasov 
type only if one uses some definite simple rules of projection. It is important that the 
same RKE can be derived also by a direct method based on the Heisenberg equation of 
motion |9|, |lO[, |16]-[^]. In the latter approach, any ambiguities in the RKE derivation 
are thus in general avoided. These features of the direct method represent a suitable 
expedient to obtain a correct generalized RKE in the frame of a relativistic analogue of 
the Kadanoff-Baym formalism. This is the main result of our study. 

Above-mentioned works |9|, |l^, |16|-|2^] are devoted to covariant generalizations of 
the Zubarev method of non-equilibrium statistical operator |47] in terms of the relativis- 
tic Wigner functions. However, this formalism presents some technical difficulties to go 
beyond the Born approximation, while the formalism introduced in this work, based on 
the real-time Green function, avoids this kind of mathematical troubles. 

In order to solve this main problem, we develop a specific relativistic modification of 
the real-time Green's functions method in the Wigner representation. This improvement 
proves to be possible since the transition to the Wigner representation in relativistic theory 
allows us to introduce a time-like direction which can be connected with the momentum 
argument of the Wigner function. As a result, the time variable determining the evolution 
of the system can be defined as a true scalar, which plays the role of the proper time in a 
given point of the eight-dimensional phase space. This invariant time is used also to define 
the chronological ordering operation in the real-time Green's functions . It leads to an 
obvious modification of this method at all stages of the formalism including the description 
of the dispersion properties of quantum field systems. Such a modification of the real-time 
Green's functions method is convenient also in order to describe the dispersion properties 
of a quantum-field system at finite temperature and density as it proves to be obviously 
relativistic invariant for all the steps of the relevant calculations. 

The outline of the paper is as follows. In Sec. 2, we discuss shortly the "direct way" 
of the derivation of RKE in terms of one-particle relativistic Wigner's functions. The 
advantages of this approach are its simplicity and clearness. This fact originates from 
the use of the Heisenberg picture from the very beginning. However, in this section, we 
do not discuss any truncation scheme and derivation of closed-form RKE. This choice 
is convenient for a comparison with the method based on the Kadanoff-Baym motion 
equations for correlation functions (Sec. 3). A covariant modification of the real-time 
Green's functions method is found as an adequate formalism for this aim. As a result, 
we obtain a generalized RKE of the Kadanoff-Baym type without using any irrelevant 
assumption, thus allowing us to use effectively the suggested method for a reliable kinetic 
description of a quite wide class of systems with inner degrees of freedom. 

In essence, the only limitation of the suggested method is a request of a polynomial 
character of interaction of systems under study. It is important that the suggested ap- 
proach results in a self-consistent, non-contradictory, and unambiguous theory of kinetic 
equations. To illustrate the possibilities of the theory, two examples are discussed in de- 
tail. Firstly, a RKE of the Vlasov type (the mean field approximation. Sec. 2) and then a 
collision integral of the BUU type (the quasiparticle approximation. Sec. 4). The standard 
Walecka model, which is a quantum field model of the relativistic nuclear matter consist- 
ing of nucleons and two types (scalar and vector) of mesons (quantum hadro dynamics) is 
used in our work as a toy model. 

In Sec. 5, we consider another generalization of relativistic kinetic theory based on 
a change of the standard averaging procedure under the equilibrium density matrix by 
an averaging using the non-equilibrium statistical operator. This step leads to a possible 
kinetic description of strongly non-equilibrium states which are often met in the relativistic 
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nuclear physics. Finally, conclusions are drawn in Sec. 6. Everywhere, we work with 
natural units, h = c = Ic. 



2 The motion equation for the Wigner function 

To illustrate our approach on a nontrivial example of a system with internal degrees of 
freedom, let us consider the kinetic description of the Fermi subsystem of a quantum-field 
system. For a consistent dynamical construction of the relativistic kinetic theory, we will 
start with the introduction of the one-particle Wigner function of the Fermi subsystem. 



Mx,p) = (2^)-^y dye-'Py < Pap{x,y) > , (1) 

where the symbol < ... >= Tr...p denotes the operation of statistical averaging with 
density matrix p in the Heisenberg representation and 

Pa/six, y) = V'/3(a;+)V'a(a;-), (2) 

V'(x) and ^(x) being the usual field operators, and x± = x i: y/2. 

In the realm of dynamical basis of the theory, we choose the Heisenberg equations of 
motion, namely, 

idf'Aix) = [A{x),P^'] , (3) 

where A{x) = A[x,'ip{x),'ip{x)] is an arbitrary local operator and P^ is the total 4- 
momentum of the system, 

pt'= [ d(j^{x\n)T^"' . (4) 



Here, day{x | n) is a vector element of an arbitrary space- like hyperplane, o"(n), with a 
time-like normal vector, n^{'n?' = 1), and T^^ = T'^^ is the energy-momentum tensor.. 

The energy and momentum conservation laws lead to the independence of the integrals 
(^) from the selection of a given hyperplane (T{n) , and we can fix the time-like direction 
by means of an external condition (a relevant example will be given below). 

In order to describe the dynamics of the system along the time-like direction, n^, and 
the space translation along the independent space-like directions on the hyperplane, (j(n), 
a boost transformation with the "velocity", n^, can be used for the motion Eqs. (|3|). 
Technically, this result is achieved using the projection of Eqs. (|3|) on the direction 
and on the hyperplane, a{n). Convolution of Eqs. ^ with the vector, n'^, leads to the 
following dynamical equation [p3|. 



i^^ = [A{x),H{r)] . (5) 

The parameter, 

T = n^x^, (6) 

plays here the role of a proper time in a new coordinate system. In Eq.(^), the derivative 
along the direction of is thus introduced, 

— - (7\ 

9r " " ^x^^ ■ 

H{t) is the scalar Hamiltonian of the system and reads as 

H{t) = n^P'' = J da{n)n^T^'''n,, (8) 
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because da^{x \ n) = n^^da{x \ n). This procedure of a formal restoration of the relativistic 
invariance in the description of the time evolution is very popular in the quantum field 
theory (see, e.g., 

The projection of Eqs. (^) on the space-like directions is carried out by the convolution 
of Eqs. (^) with the projection operator. 

The resulting equations provide a description of infinitesimal transforms of the system on 
the space-like hyperplane, fT(n). These equations are not interesting for the present aim 
of developing a kinetic theory |jl6|-|17]. 



Now, we can write the motion equation for the Wigner function (Q). Performing 
the T-differentiation of function ([l|) and using the Liouville equation in the Heisenberg 
representation, 

after substitution of the motion equation (|5|) and using the definition (0), we find: 

^^^^ = n'-^fix,?) = -i(2vr)-4| dye~'Py < [P{x,y), H{t)] > . (11) 

It is worth remarking that, for sake of convenience, the spin indices are dropped out in 
our notation. 

At this stage, we can remove the arbitrariness in the choice of the unit time-like vector, 
n^, using the specific character of the mixed Wigner's representation. Let us assume that 



the momentum vector, pf^, in Eq. (11) is a time-like one (Assumption 1). Then, we can 
define a unit vector in the direction of and identify it with the vector, n'^ (Assumption 

2) |,|ii, MM- 



= p^'/^J^ , u^ = \. (12) 

Assuming a description of non-equilibrium quantum field systems in terms of the 
Wigner functions (|l|), often one limits oneself to consider only processes on the mass shell, 

p2 = . (13) 

In this respect, introduction of Assumptiones 1 and 2 does not produce any essential re- 
striction even at the off mass shell. In other words, Assumption 1 extracts the physical 
significant domain, "P^, from the full phase space, V = {x,p). Let us remark that an anal- 
ogous restriction is usually introduced at the stage of the definition of Wigner's functions 
of the states with positive and negative energies p^]. On the other hand, Assumption 1 
can be taken off and go out into the space-like region (see the discussion at the end of 
Sec. 3). We conserve this restriction here in order to remain into the framework of the 
ordinary dynamical description. 

The effectiveness of the introduction of Assumptiones 1 and 2 is confirmed by the fact 
that substitution of the relation ([l^ ) into Eq. (11) completely removes any arbitrariness 



in this equation and leads to a "generalized RKE", 

p^d''ix)f{x,p) = -i{27rr^^Jdye-'Py < [P{x,y), H{t)] >= 

(14) 

= -ii27r)-^^ Jdye-^Py < [i/j{x+), H{t+)]'iP{x.) + i;{x+)[ip{x-), H{t_)] > , 

where we took into account the arbitrariness of the assigned hyperplane into the Hamil- 
tonian (^), H{t) = H{t^) = H{t^). Essentially, this is the first equation of the BBGKY 
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hierarchy in the Wigner representation. In the mean field approximation, Eq. (|T4|) results 
to be a non- contradictory RKE of the Vlasov type 11-0, H], H- 

In a general case, it is necessary to introduce different truncation schemes to obtain 
closed- form RKE from Eq. (|l^). The usual perturbation theory on the coupling constant 
starting from Eq. (14) was developed in Refs. 10|, [|l^]-^^. Within the framework of 
such a perturbation theory, two types of RKE with collision integrals of the second order 
have been obtained: RKE of the Bloch type (for the vertices of the "three-tails" type) and 
RKE of the Boltzmann type (for the effective vertices of the "four-tails" type). The Green's 
functions method offers an effective way to go beyond the perturbation approach. In the 
following Section, Eq.(14) will be re- formulated in terms of modified contour Green's 
functions. In particular, this will allow us to get a unambiguous generalized RKE of 
the Kadanoff-Baym type, thus opening some perspectives to go beyond the standard 
perturbation theory without breakdown of the relativistic invariance of kinetic theory. 



3 Generalized RKE in terms of modified 
real-time Green's functions 

3.1 Covariant real-time Green's functions 

The modern technique of derivation of the Kadanoff-Baym type kinetic equations is based 
on the real-time Green's functions method |]l|, ^, |3|, ^ |2^. Now, our aim is to combine 
together the " generalized RKE" (14) and a system of motion equations of the Kadanoff- 
Baym type for the real-time Green's functions. As a first step, we will modify this method 
taking into account the existence of a preferred time-like direction (Assumptiones 1 and 
2 in Sec. 2) in the definition of the Green's functions in the Wigner representation. This 
direction is fixed by relation ( [l2|) and leads to scalar time variables (see Eq. (^)), 

Ti = xfn^, n^'^u^ , (15) 

where are arguments of the two-points Green's functions. It is implied that identifica- 
tion (|l^ ) is fulfilled after a transition to the Wigner representation. The time arguments 
(|15|) are convenient for use in the definition of the ordering operation determining the mod- 
ified real-time Green's functions. Limiting ourselves to the case of a Fermi subsystem, let 
us introduce the following Green's functions in the Wigner representation, 

G^^lixp) = (2vr)-4| dy^Pyd^l{x+,x.) , (16) 

where index (k) enumerates the elements of the real-time (contour) Green's function, 

Glp{xi,X2) = -i <T''[^^{xi)i)(i{x2)]> = ^ 
Glp{xi,X2) = -i<T-[ilj^{xi)i,p{x2)]> = G++, 

G'^f^{xi,X2) = -i < llJa{xi)lpf3{x2) > = ^^2", 

G^p{xi,X2) = i < 4^p{x2)Tpa{xi) > = G^^, 

and {T"')T^ represents (anti)chronological time ordering along the r-axis: 

T-[A{xi)B{x2)\ =^ e{Ti-T2)A{xi)B{x2)-9{T2-n)B{x2)A{xi) , 
T''[A{xi)B{x2)] =^ e{T2-n)A{xi)B{x2)-9{n-T2)B{x2)A{xi) . 



(17) 



(18) 
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The operators, A{x)and B{x), are selected from the set of the field operators, ip{x) and 
ijj{x). 

The path-ordered real-time contour Green's functions are defined by the relation 



G{xi,X2) = -i <Tc[ip{xi)ip{x2)] > , 



(19) 



where Tc is the path-ordering operator on the r-contour which is chosen in the usual way 
(see Fig. T). 



n 



IrriT 



Rer 



Fig. 1. The Keldysh-Schwinger contour, C, with the upper C+ and lower C branches. 



It is worth stressing that this definition is based on the Wigner representation (|16|) 
and Assumptiones 1 and 2. This restriction leads by itself to a break-down of a simple 
analogy between the Wigner and the Fourier transforms. 

Let us return to Eq. (p^). From the definitions, (|l[) and (|l6,17), it follows that 



G^{xp) = if{xp) . 
Therefore, Eq. (0) can be rewritten in the form 

p'^a^(x)G<(xp) = (2^)-^v^ [ dye-'Py < [P{x,y),H{T)] > 



(20) 



(21) 



Let us assume now that the Hamiltonian of the fermion subsystem has a polynomial 
structure and adopt the usual notation, 



Hir) = Hoir) + F,„(r), Hinir) = Hmf{t) + F,(t) 
Here, Hq{t) describes the free evolution of the fermion field. 



Ho{t) = - I da{x\u)^p{x){--f^' (x) - m}ip{x) 



where d^{x) is the space- like derivative, 



d^{x) = A;^d,{x) 



(22) 



(23) 



(24) 



and A^'^ is the projection operator (^). In the interaction Hamiltonian, Hin{T), Eq. (22), 
the Hmf{t) part, corresponding to the mean field approximation, is explicitly separated. 
Then, Hr{T) describes the residual part of the interaction. 

It is easy to calculate the free motion part in the right-hand side of Eq. ( ^T|) for the 
Hamiltonian (p3|). For covariant transformations, it is convenient to use the following 
rule, 

'daix'\u)Six-x')yix') = ^u,ry{^), . E .(n) , (25) 



where y{x) is an arbitrary function of field operators and S{x) is the anti-commutation 



function of the fermion fields, namely S{x — x') = i['ilj{x),ip{x')]^ . Relation (25) is based 
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on a covariant generalization of the well-known property of the function S{x) (see, e.g., 

S{x) Uo=o= i7°<5^'na;) , (26) 
(a simple proof of Eq. (25) is given in Appendix A). This kind of calculation for Eq. ( pl| ) 
leads to the following result: 



p^^df^{x)G^{xp) + ^Ip^f^ , d^{x)G^ (xp)] + im[p,G^{xp)] = 

= {27r)-*^ J dye-'Py < mx+), HUr+M{x-) + i^{x+)mx^), H,^{t_)] > , 



(27) 



where p = p^j^. 

The following step consists in the search of a connection between the right-hand side 
of Eq. (|27| ) and the mass operators, T,^''\xi,X2), defined on the r-contour. To this aim, 
let us write the motion equations for the real-time Green's functions, 



[ijd{xi) -m] G{xi,X2) = Sc{xi,X2) + J^d'^x' T.{xi,x') G{x',X2) 
G{xi,X2) [-i"f d {x2)-m] = 5c{xi,X2)+ J^d'^x' G{xi,x') S(x',X2), 



(28) 



The integration in the right-hand side of Eqs. (^) is performed according to the rules: 



d^x... 



and 



p roc rex 

/ dr... = / dr... \^+ - / 

J C Jt() J to 



da{x\u) / dr.... 



dr.. 



(29) 



(30) 



where c"'"(c ) is the upper (lower) branch of the r-contour (see Fig. 1). Finally, the delta 
function on the r-contour is defined as 



Scixi,X2) = Saixi,X2)5ciTl -T2) , 



where 



^cin -T2) = < 



5{ti — T2) if ri G and T2 £ c'^ , 
-6{ti — T2) if Ti £ c~ and T2 £ c~ , 
otherwise . 



(31) 



(32) 



5a{x) is the d=3 delta function of the space arguments on the hyperplane a{x\u). 

For purpose of comparison with Eq. (p7[), let us write the motion equations for the 



correlation function, G^(xi,X2). From Eqs. (^8|) , the Kadanoff-Baym type equations 
follow 



[i7a(xi)-m-SMF(xi)] G<ixi,X2) = I dx'{J:<{xi,x') G^(x', X2)+S«(xi, x') G<(x',X2)} 

(33) 



G<(xi,X2) Ht d {x2)-m- ^mf{x2)] = / dx'{ G<{xi,x') S^(x', X2)+G^(xi, x') S<(x',X2)} . 

(34) 

Here, the mean field part of the mass operator was extracted according to the following 
prescription: 



S(xi,X2) = T,MFixi)5{xi - X2) + 

+ 0(ri-r2)S>(xi,X2)+e(r2-ri)S<(xi,X2) . 



(35) 
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Moreover, the retarded and advanced Green's functions were introduced, 



G^M(X1,X2) = ±^[±(ri-T2)]{G>(xi,X2)-G<(xi,X2)} = 

(36) 

= G%X^,X2) - G</>(X1, X2) = G>/<(X1, X2) - X2) . 

In a similar way, the retarded and advanced components of the mass operator are de- 
fined. A mutual exchange of the symbols > and < transforms Eqs. ( p3D and (^) into 
the corresponding equations for the correlations functions, G'' {xi,X2)- By this way, the 
Kadanoff-Baym equations (|3^ and (|3j) have the usual form however they are defined 
using covariant modified constituents (Green's functions and mass operators). However, 
these small corrections are essential for the search of the agreement with RKE (p7| ) . 

The transition to the mixed Wigner's representation for the motion Eqs. (^) and 
( |3^ is usually done with a limitation to the lowest orders of the gradient decomposition 
(the first order into the drift part of RKE and order zero for the collision part) p9[ |, 
approximation that is justified only for rather slow kinetic processes. This allows us to 
write the motion equations in the simplest local form. This restriction casts some doubt 
upon the description of non-equilibrium processes in nuclear matter at extreme conditions 
(see, e.g., pT] , [28[|). However, for simplicity's sake, we keep here this approximation only. 
Then, the motion Eqs. (|3^ ) and (|3^ in the Wigner representation are transformed into 
the following equations: 



o^^x): 



\p+l^d{x) - m - Smf(x) + ^d^{x)^MF{x)d''{p)]G<{xp) 



^<{x)G'^{xp) + T."-{x)G<{xp) , (37) 



G<{xv)[v - ^7 5 (x) - m - Smf(x) " ^ ^ {p)d^'{x)T.MF{x)] 



G<(xp)S^(xp) + G^(xp)S<(xp) . (38) 



We have neglected gradient terms in the right-hand sides of these equations [29]. 



3.2 Generalized RKE 

Let us fulfill now the concordance of the " generalized RKE" (|l4|) and the system of equa- 
tions of the Kadanoff-Baym type ( pT] ) and ( p8| ) . Obviously, for this aim, it is necessary to 
find such a combination of these equations which provides a coincidence with the left-hand 
side of Eq. ( [27| ) describing the free evolution of the fermion subsystem. This assumption 
will permit us to identify (in the corresponding orders of the gradient expansion) also the 
other parts of Eq. (|^ and Eqs. (|^), (||). 



In order to perform such a construction, it is sufficient to multiply Eq. (pq) by ip from 
the right and Eq. (^) from the left and to subtract the second result from the first one, 
i.e., it is necessary to consider the following algebraic combination: 



Eq.{2,S)ip-ipEq.(^ . (39) 
Then, we get the following equation 

p^'d^{x)G<{xp) + ]^[pipi,d^{x)G<{xp)] + im[p,G<{xp)] - 

-i{G<{xp)T.MF{x)p-pT.MF{x)G<{xp)} + (40) 
+\{d^{p)G<{xp)d^'{x)T.MF{x)p + pd^{x)T.MF{x)d^'{p)G<{xp)} = -S{xp) , 
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where S{xp) is the cohision integral, 

S{xp) = ip{^<{xp)G^{xp) + ^^{xp)G<{xp)} - 

-i{G< (xp)S^(xp) + G^{xp)^< {xp)}p . (41) 

Now, from a comparison of Eqs. (|2^) and (|^), it follows an equivalence of the left and 
right parts of those relations: 

(27r)-^y^ J dye-'Py < [i;{x+), HMF{r+mx_) + i;{x+)mx_), HMF{r-)] >= 

= i{G<{xp)^MF{x)p - p^mf{x)G<{xp)} - (42) 



-{d^{p)G<ixp)d^'ix)^MFix)p + pd^ix)^MFix)d^'{p)G<ixp)} 



Sixp) = -(27r)-Yp'y dye-'Py < [V;(x+), F,,(r+)]^(x_) + ^(x+)[V'(x_), /f,,(T_)] > . 

(43) 

Let us remind that these equalities are correct only in the lowest orders of the gradient 
expansions. 

Both parts of the equality ( |4^ ) describe the force contribution (with the inverse signs) 
of the mean fields into the convective part of RKE (^^ or (^0|). Let us illustrate this 
aspect by means of an actual model. We thus refer to the standard Walecka model, where 
the Hamiltonian for the interaction of nucleons with massive vector and scalar meson 
fields reads as [30|: 

Hin{T)= da{x\u)i){x){gyUJ^'{x)-f^- gs(l){x)}'4){x) . (44) 



In the mean field approximation, the relevant Hamiltonian, HMFiT), can be deduced from 
this one by replacing the meson fields operators with their mean values, (p ^< (j) >, and 
uj^ — >< io^^ >. Then, the corresponding contribution to the mass operator is equal to 

^mf{x) =gv< w^(x) >jf,-gs < (t>{x) > . (45) 

In the first order of the gradient expansion, both RKE (^) and (^0|) lead to identical 
results in the frame of the mean field approximation: 

P''d^{x)G<{xP) + ^d''{x)M{x){P,d^{P)G<{xP)} + 
+g,P''F^4x)d''{P)G<{xP) + ^[P-f'',d^{x)G<{xP)] + (46) 
+iM{x)[P,G<{xP)]-y,F^,{x)[P-f'',d''{P)G<{xP)] = -S{xP) , 

where M(x) = m — gg < 4>{x) > is the effective nucleon mass in the mean field approx- 
imation, F^^{x) = < u)'^ > —d'^ < uj^ > and P^^ = p^ — g^ < uj^^{x) > is the kinetic 
momentum. 

RKE ( pBD converts in RKE of the Vlasov type when conditions of the collision effects 
are neglected {S{xP) = 0). The Vlasov RKE in this form was obtained and analyzed 
for the first time in Refs. [|l6|-|17]. These equations are rather complicated because the 
spin and meson degrees of freedom as well as states with positive and negative energies 
are taken into account. On the basis of this RKE, a series of more simple, particular 
cases can be analyzed. The simplest situation corresponds to the case of a spin saturated 
system without antinucleon states, where all the spin-dependent effects in the RKE can 
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be neglected (the spin effects were investigated separately in Ref. [0). Using the Clifford 
decomposition of the Wigner function, we get - in the quasi-classical limit - the well- 
known RKE of the Vlasov type in the quantum hadrodynamics, 

P^'d^{x)f^ixP) + Md''{x)Md^{P)f{xP) + g,P^F>'''d,{P)f{xP) = , (47) 

where f^{xP) corresponds to the scalar part of the Clifford decomposition of function 
G<{xP) (see Eq. (|o|)). It is important that the basic RKE (H) leads to a correct RKE 
( p7| ) of the Vlasov type for a spin saturated system without any additional assumption, 
like that of the rescaling procedure type (see the detailed discussion in fl^). 

Technical details of derivation of Eqs.(^)-(|47|) are given in Refs. It is not 

difficult to obtain analogous RKE also without use of the gradient expansion. Such a 
RKE has a non-local character in the momentum space. 

Calculation of collision integrals can be performed using two alternative schemes based 
either on formula (HTl) or (03). In the latter case, a standard perturbation theory was 



developed in Refs. ||9|, |l^, [l^, 1^, 2C]. Within the framework of the Walecka model ([44[), 
this approach allows us to derive collision integrals of the Bloch type (when the meson 
subsystem is considered on equal terms with the nucleon one) or of the Boltzmann type 
(if the role of mesons is confined only to the formation of the effective nucleon-nucleon 
interaction) . 



Let us return to the representation of the collision integral (41) based on the method 
of modified real-time Green's functions. Our subsequent transformations are based on the 
analogy with the nonrelativistic theory |]2^ , p9| . From the definition (|36|), it follows that 
the Green's functions and mass operator satisfy the identities: 

G^ixp) - G^ixp) = G>{xp) - G<{xp) , ^^g^ 
Ti^{xp) — T,'^{xp) = T,''{xp) — T,'^{xp) . 

Using these relations, the collision integral (^) can be transformed into the form 



(49) 



S{xp) = ^{p[J:>{xp)G<{xp)-J:<ixp)G>{xp)] + 

+ [G<{xp)J:>{xp) - G>{xp)^<{xp)]p} + S°ff{xp) , 

where ^ {xp) is the off-mass-shell part of the collision integral, 

S°ff{xp) = ^{pi:<{xp)[G^{xp) + G^{xp)]+p[T.^{xp) + i:^{xp)]G<{xP)- ^^^^ 
- G<{xp)[T.^{xp) + T.^{xp)]p - [G^{xp) + G'^{xp)]T.<{xp)p} . 

Of course, the off-shell effects enter in the whole collision operator (49), but function 
(50) refers to them exclusively for the present aims. The separation into two parts, (49) 
and (50), is ascertained by means of spectral properties of the Green's functions and 
the mass operator. In our formalism, the corresponding spectral decompositions can 
be fulfilled in a covariant form (see Appendix B). This circumstance is a very desirable 
property of the relativistic kinetic theory which is intended for a description of macroscopic 
motions of a medium. Let us remark that the off-shell processes play the most important 
role in the relativistic nuclear physics of intermediate p9|] a nd high |^0| energies. 

The generalized RKE ( ^0|) with the collision integral (|49|), (|50| ) are the main results of 
this Section. It is worth remarking that a few other versions of generalized RKE have been 
previously obtained in the literature (see, e.g., and |@). We have in mind the cases of 
derivation of RKE for the Fermi subsystem based either on the Dirac type equations of 
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motion for the field operator or on the equations of the Kadanoff-Baym type (p^), (^4|). 
However, it is clear that this procedure contains some elements of an indetermination 
and arbitrariness on the stage of the construction of the motion equation for the Wigner 
function. In considered our case, the selection rule (^) has a reliable dynamical basis 
(see Sec. 2) which was tested in Refs. [17|-|22]. Of course, the generalized RKE (40) can 



be obtained also by a direct way on the basis of Eq. (P7|), without using the rule (|39[) , 
with a help of the Dyson equation and definitions of the type of Eqs. ([T7|). 

RKE (po|), (|49|), and ( ^0|) have a more complicated matrix structure. This character- 
istic behavior arises by taking into account inner degrees of freedom which are intrinsic 
for the considered Fermi subsystem and, in particular, states with positive and negative 
energies. Indeed, the presence of 7^ matrices before the time derivative in equations of 
the Kadanoff-Baym type ( p3| ) and (p^) ensures the "direct" time movement for states of 
the spinor fields with positive energies and the "inverse" time direction for states with 
negative energies. But for the kinetic theory, it is characteristic the existence of a single 
time arrow. The rule (|39| ) just ensures a covariant way of projection of Eqs. (^) and (^^ 
on such a single time direction. 

Finally, let us return once more to the role of the assumption > (Assumption 
1). A consecutive realization of this limitation leads to appreciable complication of the 
considered integrals. Similar restriction is absent in the traditional approach where an 
integration is carried out in the whole momentum space. In order to get traditional results 
in the frame of our formalism, it is necessary to implement the extension of obtained RKE 
into space-like region. There are three aspects of such a continuation: 

1) into the integrals with respect to the momentum space, 

2) into the time differential operations, 

3) into the time integral operations. 

The extension of the momentum integrals into the region, p'^ < 0, is trivial in the 
quasiparticle approximation (Sec. 4). 

As for the derivation of the Vlasov RKE, the considered continuation does not lead to 
some changes of differential operations with respect to the time (e.g., to a change d/dr 
by id/dr). Indeed, in the mean field approximation, the calculation of the force part of 
the Vlasov RKE (the right-hand side in Eq.(p7|)) results in the cancellation of once 



rule (25) has been used. In the other words, RKE conserves its physical content. 

Finally, in the Markovian approximation, integral time operations do not lead to any 
modification of results at the continuation into the region j?^ < 0. 

The possibility of continuation beyond the above-mentioned cases needs further re- 
searches. 



4 Quasiparticle approximation 

In this Section, we derive a covariant generalization of the quasiparticle approximation 
(QPA), widely used in the kinetic theory, which allows us to leave, by the simplest way, 
the framework of the usual perturbation theory. After that, we check possibilities of the 
suggested method on an example of derivation of collision integral of the Boltzmann- 
Uehling-Uhlenbeck (BUU) type for the standard Walecka model of relativistic nuclear 
matter consisting of nucleons, scalar and vector mesons with the interaction Hamiltonian 
(^^. Finally, we discuss briefly a feasible generalization of the QPA in the model (p^). 

Let us introduce a non-equilibrium spectral function in the Wigner representation (see 
the first formula in Eq.(48)) 

a{xp) = i[G^{xp) - G^{xp)\ = i[G>{xp) - G<{xp)] , (51) 
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This function satisfies the following sum rule 

(27r)^y" dEa{x;p_L,E) = ju , (52) 
which is a simple consequence of the single-time anti-commutation relations for the field 



operators and of the rule (25). In Eq. (|5^), the integration is carried out over the 
longitudinal energy, E = pn; pj_ = /^^'^pu is the transverse momentum (A^*^ is the 
projection operator (^)). 

G^{xp) and G'^{xp), can be represented as: 

G<{xp) = ia{xp)T{xp) , ^^^^ 
G^{xp) = —ia{xp)[l — T{xp)] , 

where J-{xp) is a unknown function. 

To determine the non-equilibrium spectral function, we consider relation ( [5l| ) and the 
corresponding motion equations for the retarded and advanced Green's functions. We 
write at once these equations in the minimal order of the gradient expansion, namely, 

\p + ^-fd{x) -m- ^mf{x) + i9^(:E)SMF(x)a^(p)]G^/^(xp) = 

= l + S^/^(xp) G^/^{xp), 

G^M(xp)[p - §7 a (a;) - m - ^mf{x) - | dAp)9nx)^MF{x)] = 

= l + G^/^(xp) S^/^(xp). 



(54) 



Since the gradient expansion orders in the drift and collisions parts of RKE (|40| ) are 
fixed, we can restrict ourselves to order zero of this decomposition in the estimation of 
correlation functions (^), so that 

ip-m-J^MFix) - S^/^(a;p)] G^/^(xp) = 1 . (55) 

In our toy model, T,mf{x) is defined by Eq. (|45|) and hence, 

[P - M(x) - S^/^ (xp)] G^l^ {xp) = 1 . (56) 

Let us assume now — as a working approximation — that the mass operator, E(xp), has 
the same matrix structure as the mean field part (|45| ) (in the Walecka model (^^, this 
implies a possibility to neglect the tensor part of Clifford's decomposition of the mass 
operator) , 

E(xp)^S,(2;p)+7^E^(xp) (57) 

(a more general form is given in Ref . ) . In this expression as well as the following ones, 
the marks R/A are omitted for simplicity's sake. Eq. ( [5^ ) has then a matrix structure 
like the Dirac equation of noninteracting fields (quasiparticle structure) 

[V{xp) - M{xp)] G^l^ixp) = 1 , (58) 

where 

'P^l{xp) = P^i{xp) - D^(xp) , (59) 
M{xp) = M{xp) - ^s{xp) ■ (60) 
Therefore, Eq. (^) can be easily solved: 
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Then, from Eqs. ( |5lD and (|6l|), we get 

a{xV) = 2tt{V + M)5[Re{V^ - M^)]{9{E) - 9{-E)} . (62) 
This result is correct, provided that the quasiparticle excitations are weakly damped, i.e., 

\Re{V^ - M^)\ » I Jm(p2 _ _^2)| _ (gg^ 

It is easy to verify that the spectral function (p^) satisfies the sum rule (^). 

The substitution of the non-equilibrium spectral function ( |62|) into Eqs. ([5^ ) leads to 
the following expressions for the correlation functions. 



G<{xr) = 2Tri{V + M)6{V'^ - M'^)J'{xr') , 

G> {xV) = -2m{V + M)5{V'^ - M''^)[l - J'ixV')] , 



(64) 



where V' = ReV and M' = ReM. 

As an illustration, we present a derivation of a collision integral of the BUU type 
in the Walecka model (^^. This example is well known in the literature, therefore it 
is convenient for testing our method. We confined ourselves to the Born approximation 
only in the calculation of the self energy parts, E^/*^, in the collision integral (|49|). The 
corresponding diagrams are shown in Fig. 2. 









Fig. 2. Born diagrams for the self-energy part, — iS^ = — The solid, dashed 
and wavy lines denote the propagators of nucleons, scalar and vector mesons, respectively. 



In a general case, this collision integral has a rather complicated form since we take into 
account simultaneously the spin and meson degrees of freedom and states with positive and 
negative energies. The simplest situation corresponds to a spin saturated nuclear matter 
without antinuclear component. In order to get the corresponding collision integral, we 
have to perform a transform from the spinor representation to the spin one for states 
with positive energies ||2^ and to take into account only the diagonal part of the Wigner 



function in the spin representation. As a result, for function J^{xV) in Eq. (64), we obtain 
(for simplicity, we omit here and below all primes in functions M' and V') 



E 

r,s=l,2 



ixV) 



2po 



<^{xV) , 



(65) 
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where ^srixV) is the Wigner function of states with positive energy in the spin repre- 
sentation. For a spin-saturated nucleon subsystem we have — ^^rs, where ^{xV) is 
the scalar Wigner function. In Eq.(^), we used the completeness relation for the Dirac 
spinor [^] 

v-^ r + M 

2^ uu - 



Substitution of Eq. ( |65[) into the first formula of Eqs. (64) leads to the relation 

Q<{+) (^p) = iT,'^l±^s(v° - VM^ + V^^ F{xV) , (66) 
where a new auxiliary function is introduced: 

^(^^) = ^^(^^) • (67) 

This expression differs from the true Wigner function, ^(xV), by the scaling factor, A4 /V^ 
[42|. Introduction of such a rescaled Wigner function (^) is convenient in order to write 
the collision integral. 

A tedious but straightforward evaluation of the on-mass-shell part of the collision 
integral (E^) leads to the following result: 



{F{xV)F{xV2)F{xV3)F{xVi) - F{xV)F{xV2)F{xV3)F{xVi)} , (68) 

where F = \ — F is the Pauli blocking factor. Here, the transition rates are defined as 
follows: 

W = Wss + W,, + Wsv ■ (69) 

The partial transition rates are given by exchange of the corresponding pairs of the scalar 
or vector mesons: 

HVV:,){V2Va)\ - Ds{V - r3)Ds{r - V^)[M^ + 
+M'^{PV2 + rVs + VVa + 7^2^3 + 'P2V4 + VsVa) + 
+{'P2Vi){VV3) + {V2V3){VVi) - {VV2){V3Va)]} , 
W,,{VV2V3Va) = 8gt{ADl{r-r3)[2M^-M\V2V4 + rr3)+ (70) 

+{'PT2)ir3r4) + {vva){V2V3)] - 

+PP3 + VV2 + V3VA + ^27^4) - 1{VV2){V3V4)\} , 

Ws^imVsn) = 8glgl{-ADs{V - V3)D,{V - V3)M'[VV2 + 

+VVi + V2V3 + V3V4] + Ds{V - V3)D,{V - Va)[2M^ + 
+M'^{-VV2 + 2VV3 - VVa + 2V2VA - V2V3 - 
-V^Va) + 2{VV3){V2Va)]} . 

Finally, Ds^vi'P) are the Fourier transforms of the mesons Green's functions, Ds^y{V) = 
(P^ — rri^ , and and rriy are masses of the scalar and vector mesons, respectively. 
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It is worth noticing that similar results have been recently obtained in Ref. [Q] within 
the framework of non-covariant real-time Green's functions method. A comparison of our 
results with those of Ref. Q shows an overall good agreement (up to a replacement of the 
fermionic Green function calculated in the mean-field approximation by that obtained 
in the more general approximation (64)). However, it needs to keep in mind that the 
standard Walecka model is too primitive for a proper description of real properties of 
nuclear matter at intermediate energies. Already at an equilibrium state, it leads to 
an equation of state stiffer than the expected one and, at moderately high density and 
temperature, the effective masses of nucleons become very small or even negative |33]. 
At present, two directions of the Walecka model improvement exist: either an increase 
of the number of its constituents [^] or a transition to a nonlinear generalization of the 
model [ p3| , |35[ . The first way leads to technical complications of the kinetic theory (by 
virtue of an essential increase of the number of constituents for real system), the second 
one remains still unexplored for a description of non-equilibrium states. 



5 Connection with the non-equilibrium statistical 
operator method 

As previously outlined, in the definition of Wigner's (||) and Green's functions (0), the 
statistical averaging was performed by means of the equilibrium density matrix, p, in the 
Heisenberg picture corresponding ordinarily to a system state at an infinite past. This fact 
limits the applicability of the theory to the case of weakly non-equilibrium states which 
can be described defining only slowly varying thermodynamic functions of the system 
such as temperature, chemical potential, etc. Now, let us discuss briefly a possibility of 
a generalization of the formalism suggested in Sec. 3 which is based on the change of 
the statistical averaging procedure, using the non-equilibrium statistical operator, p{t). It 
is expected that such a generalization will allow us to describe strongly non-equilibrium 
states of a system p4| , p9| , |3^ . The non-equilibrium statistical operator method was 
adapted to fulfill requirements of the relativistic kinetic theory in our previous works 
p!o| , p!^ , [r9| . Here, we discuss only a problem where this method is combined with the 
covariant formalism of real-time Green's functions method (Sec. 3). 

Let us write the motion equation for the non-equilibrium statistical operator in a 



differential form |1C] 



dp{T)/dT = -e{p{T) - p,{t)} , (71) 

and, in an integral form (e is an infinitely small value, e > 0, and will approach to zero 
after the execution of the thermodynamic limiting transition), 

T 

p{r) = p,{T)+i J dT'e<^'-^^ {[p,{T'),H{r')]+idp,{T')/dT'] . (72) 

— oo 

Here, the time, r, is defined by Eq. (^) with a subsequent introduction of Assumptiones 
1 and 2 and relation ([l2|). An infinitely small source in the right-hand side of Eq. ( |7T|) 
is introduced to break the symmetry of this equation with respect to the time refiection 
(in the Wigner sense). The form of Eq. ( [Til) allows us to select the retarded solution of 
the Liouville equation on the basis of an analogy with the formal scattering theory. This 
implies that after completing our calculations the thermodynamical limit should be taken 
and, then, the limit e — > assumed. 



The quasi-equilibrium statistical operator, Pqir), in Eq. (|7lD is the asymptotic form of 



the non-equilibrium statistical operator, p{t), at the kinetic stage of evolution, when r 
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— oo. An explicit form of operator Pqir) can be derived from the principle of the maximum 
of the information entropy under a supplementary condition of the given averaging value 
of < Papixy) >qT, where < ... >qr= Tr...pq{T) and the operator, Pa/3{xy), is defined by 
Eq. (HI). For our task, it is convenient to associate the proper time, r, with the slow 
variable according to the expression (^. 
As a result, we have |10| 

Pg{r) = exp{-S{r)} , (73) 
where S{t) is the entropy operator of the system at the kinetic stage, 

S{t) = <^{t) + J da{x\u) J dS Pap{x,y) F^p{x,y), (74) 

and <&(t) is the normalizing functional defined by condition Tr Pq{T) = 1. The Lagrange 
factor, Fai3{x,y), is determined by the self-consistency condition 

< Papixy) >T = < Pa/sixy) >gr ■ (75) 

In the left-hand side of this equality, the averaging is performed by means of the non- 
equilibrium statistical operator p{t), i.e., < ... >t= Tr...p{T). Relation ( |75[) shows a full 
equivalence of both statistical operators, p{t) and Pqir), at the kinetic stage of evolution. 
Hence, both variants of the Wigner function. 



fap{x,p) = (27r)-^y dye-'Py < Pap{x,y) >r= (2^)-"y dy e'^Py < Pap{x,y) , 

(76) 

can be used for a kinetic description of both quasi-equilibrium and non-equilibrium states. 

The Green's functions, (|T^)-([T9|) and (^6|), also permit an analogous generalization for 
the case of strong non-equilibrium states. 



Ginon-eq){xp) = -i(27r)-^ / dyc'^^ < r,[V'(x+)V(x_)] . (77) 



The form of the Kadanoff-Baym Eqs. (p3[), (p^) remains also the same, taking into ac- 
count that the e-terms generated by the right-hand of Eq. ( |72[ ) vanish after a realization 
of the thermodynamical limit. However, a restriction to the lowest terms of the gradi- 
ent decomposition is incompatible with a description of strong non-equilibrium states. 
Therefore, the corresponding results of Sec. 4 require a further generalization. This field 
of relativistic kinetic theory is still scarcely investigated. 

Let us remark that Eq. ( [72| ) can be rewritten in a form convenient for the construction 
of an iteration scheme under the energy of particle-particle interactions. In full analogy 
with the non-relativistic case |2^], we obtain |^, |l^ 

p{T)=p,iT)+i j dT'e<-'-^^{[p,{T'),HUr')] + 

— oo 

+ /da(x»/dV < [Paf^{x',y'),H,^{r')] .^J^fS].^ } ■ (78) 



As shows our experience, the second addend in the right-hand side of Eq. ( |78D removes 
usually all disconnected diagrams of a given order from the set of diagrams which are 
supplied by the first addend. 

Some problems of this approach arise on the level of the non-equilibrium thermody- 
namical Wick's theorem. However, there is a set of methods which permits to overcome 
these difficulties |4^ . c 
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6 Conclusions 



In the present work, we have demonstrated an interesting possibility for a covariant gener- 
ahzation of the real-time Green's functions method which opens new ways after a transi- 
tion to the Wigner representation. In fact, this generalization is based on the admissibility 
of introduction in every point of the Minkowski space a unit time-like vector constructed 
with the momentum vector of the corresponding point of the phase space. Therefore, 
for a definition of the preferred time-like direction in the Minkowski space, we used a 
possibility contained in the Wigner function itself. The relevant Assumptiones (Sec. 2) 
allow us to write, starting from the Heisenberg representation, a motion equation in a 
covariant form for the one-particle Wigner function (Sec. 2) and, after that, to obtain 
a generalized RKE in the frame of the Kadanoff-Baym technique (Sec. 3). Such a way 
is found especially effective for a kinetic description of subsystems with inner degrees of 
freedom and allows us to eliminate some ambiguities of relativistic kinetic theory based 
on the Kadanoff-Baym equations. 

As a test of the suggested approach, the Fermi sector of the well known Walecka model 
of the relativistic nuclear matter was chosen (Sec. 3 and 4). For this model, RKE of the 
Vlasov type (Sec. 3) and of the BUU type (Sec. 4) were obtained. In the latter case, the 
collision integral was derived only for a spin saturated system without antinucleon states. 
The calculations were performed within the framework of the quasiparticle approximation 
(QPA). We discuss also a possibility of an extension of this approximation. An agreement 
with a set of well-known simple results in the literature denotes a noncontradictory char- 
acter of the proposed generalizations. It is safe to assume that methods suggested here 
provide also a correct description of spin degrees of freedom, states with negative ener- 
gies (e.g., for the annihilation channel in the theory of non-equilibrium electron-positron 
plasma), and so on. 

Finally, the covariant modification of the real-time Green's functions method suggested 
in Sec. 3 was generalized to a kinetic description of strongly non-equilibrium states (Sec. 
5). For this aim, we suggest to use the non-equilibrium statistical operator method. 

The non-equilibrium relativistic nuclear matter represents a good example of very 
complicated object for theoretical and experimental investigations. Of course, the simplest 
toy model and our present results discussed in Sec. 4 are only a rough approximation of 
the reality. At present, a lot of actual problems still remain in the relativistic kinetic 
theory, e.g., going beyond the quasiparticle approximation, taking into account nonlocal 
(and, in particular, non-Markovian) effects in kinetic processes, cluster decomposition, 
etc. (see, e.g., ||l], ^ 37, |3^, 43, 44, ^). An important feature of the examined 



approach is a lack of some dynamical constraints into a corresponding quantum field theory 
(these constraints are introduced in order to eliminate " unnecessary" , non-physical degrees 
of freedom). The dynamical introduction of these constraints into the Kadanoff-Baym 
formalism is a rather non-trivial problem. A necessity to consider such constraints creates 
serious problems for the kinetic theory (e.g., in order to describe the A-isobar subsystem 
of nuclear matter). From our point of view, the approach suggested in the present work 
can serve as a reliable basis for such kind of further researches in the relativistic kinetic 
theory. 
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Appendices 
A Proof of Eq. (25) 

Let us derive a proof of the Eq. ( [25| ) using a boost transformation. 

I d^z'Siz-z',0)fiz') = ij^fiz) , (Al) 

for an arbitrary function of the field operator f{z). The integration here is fulfilled over 
hyperplane of the constant time, z^ = z'^. Let us perform a transition to an arbitrary 
frame of reference using a transformation from the homogeneous Lorentz group 

= AtM^^ x'" =AtH^'^ {A2) 

The velocity vector, u^, defines simultaneously the orientation of the hyperplane, a{u). 
The transformation of the commutator function, S{z — z',0), in Eq. (Al) is obtained 
using the corresponding unitary operator, U{A) 



S{z - z', 0) = S[A{x - x')] = U{A)S{x - x')U~\A) , (A3) 

where x,x' G a{u). Similarly, we have: 

fiz) = UiA)fix)U-\A) . (AA) 

Let us take into account also the relation 

U-\A)-f^'UiA) = At(n)7^ . (A5) 

This equality is a consequence of the invariance under the Lorentz transformation (A2) 
of the scalar product of ^p{x) and 0(x) p9[|. 



{ip,(l))a = J da^{x\u)'il:{x)-if,(t){x) . 

Finally, it is necessary to consider the equalities A'q(u) = and A']^(n) = u^{k = 1, 2, 3) 
Formulas (A2)-(A3) allow us to realize the transformation from Eq. (Al) to Eq. ( p5|) 



Another proof of Eq. (25) can be found in Ref. jlO|. 



B Spectral properties of the modified Green's 
functions 

The spectral representations of the covariant contour Green's functions introduced in Sec. 
3 can be also considered in a covariant form (we follow the work Q were a similar approach 
was fulfilled in the framework of a non-covariant formalism). The main feature of these 
representations is the selection of the time arguments in the form (15). As an example, let 
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us consider the spectral decomposition of the retarded and advanced Green's functions. 
According to the definitions, (f6[) and (|36|), we have (see the last remark in Sec. 3.2) 



G'"^{xp)=±{2^)-^ j dye'^y{G>{x+,x^)-G<{x^,x^)}e[±y,u^] ., 
where is the unit vector (^). From here, it follows 

CO 

G>{xp') -G<{xp') 



G^/A(xp) = ±— I dV / dr/5W(p^-p''^ + n^r/)- 
2m J J 



±r] — ie 

— OD 

or 

G^/^ixp) = ±^[G>{xp) - G<{xp)]+ G [xp) , (Bl) 

where 

oo 

G (xp) = ^-j A' j dv 6^^^ (p^ - p'" + u^ri) [G> {xp') - G< {xp')]V {^^ . {B2) 

— oo 

Let us transform the last relation using the decomposition of an arbitrary vector, a^, on 
the time-line direction, (12)) and orthogonal to it space-like component, 

a^' = uf'iau) + Ai^^a" = af + , 

where the projection operator, A^'^, is defined by Eq. (^). Then, we can write the 
following representation 

5(4)(p/. - p't" + u'^r]) = 6[n^{p'' - p'f" + u'^T])] (5^(p^ - p^) , {B3) 

where function 5o-(p^) is defined on the space-like hyperplane a{u) . Using Eq. (B3), we 
can rewrite Eq. (B2) in the final form: 

oo - ^ 

G (xp) = j dojV { ^ \ {G>{x-p^,oj) - G<{x;p^,u;)} . {BA) 

— CO ^ ' 

The denominator, {uj — Vp^), in Eq. (B4) shows clearly that function G {xp) is just the 
off-mass-shell part of the Green's functions, G^/^{xp). As far as we have 

G^{xp) + G'^{xp) = 2G{xp) , 

and an analogous relation for the mass operator, S^/"^, the last term in Eq. (^) represents 
indeed the off- mass-shell part of the collision integral (^8|). 
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